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in the presence of a strong brane-bulk coupling is vastly different from the 
usual case of a weak coupling. In particular, some active neutrinos could be 
absorbed by the bulk when they oscillate from one kind to another, a signature 
which can be taken as the presence of an extra dimension. In a very large class 
of models which we shall discuss, the amount of absorption for all neutrino 

X' 

$—1 ' oscillations is controlled by a single parameter, a property which distinguishes 

extra dimensions from other mechanisms for losing neutrino fluxes. 



I. INTRODUCTION 

In the braneworld scenario [1], only gravity and sterile neutrinos can move off our 3- 
brane into the bulk, so it is only these two kinds of objects that can be used to probe the 
extra dimensions. If the size of the extra dimensions is considerably smaller than 0.1 mm, 
or there is at most one extra dimension with that size, gravity cannot be used to reveal 
the presence of an extra dimension, so it is up to the neutrinos. Assuming the coupling 
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between the neutrinos in the brane and the bulk to be large, we shall argue that it may be 
possible to detect relatively small extra dimensions using precise neutrino oscillation data 
in the (distant) future. 

The strength of this coupling is of course unknown, but in our minds it is likely to be 
large. It is usually assumed that the coupling between the active brane neutrinos and the 
Kaluza-Klein (KK) neutrinos of the bulk is of the Dirac-mass type. If this Dirac mass is 
comparable to the charged fermion masses, then the brane-bulk coupling is indeed strong 
because the measure of its strength is the ratio of the Dirac mass to the Majorana mass of 
the neutrinos, a very large ratio. 

This coupling is taken to be small in most of the discussions in the literature [2]. In 
that case the result follows from perturbative calculations. Some, such as Barbieri et al, and 
Cosme et al [2], also considered the effect at large couplings via numerical simulations. In 
Ref . [3-5] , we studied the strong coupling limit of a minimal model (MM) analytically. This 
is a five-dimensional model with the smallest number of parameters which can accommodate 
the present experimental results on neutrino oscillations. The details and the consequences, 
which will be explained more fully below and in the next section, turn out to be very different 
from what one can expect by extrapolation from the weak coupling limit. In particular, to 
get three mass eigenstates on the brane in the strong coupling limit, we have to start out 
with four brane neutrinos when the coupling is turned off. The results are also distinct. A 
small solar to atmospheric mass-gap ratio naturally leads to a small reactor angle # 13 , as 
observed. There is also an absorption of the active neutrino fluxes by the bulk. We shall 
argue below that this absorption, if experimentally observed and if it satisfies the properties 
to be described later, is a good indication of the presence of extra dimensions. 

The purpose of this paper is to show that these distinct predictions are generic, whose 
validity goes way beyond the MM. We shall produce a large class of models, including those 
with a non-zero bulk masses, and those with non-trivial bulk-bulk interactions, as well as 
some models with more bulk neutrinos and/or in more extra dimensions, which yields the 
same essential properties as the MM, and gives rise to the same predictions. This finding 
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then allows us to accept the consequences as a generic feature of the extra dimensions, rather 
than some artifacts coming from the specific assumptions of the MM. 

We shall now define the MM and summarize its properties in the strong-coupling limit. 
For more details please consult the next section. The MM is a model in five spacetime 
dimensions, with a single massless neutrino in the bulk. Its KK states have an integral mass 
spectrum, in an energy unit which is inversely proportional to the size R of the extra dimen- 
sion. We shall adopt this unit throughout, so that all masses are expressed as dimensionless 
parameters. There are four left-handed neutrinos on the brane, the usual three, plus a ster- 
ile neutrino. They are assumed not to interact directly among themselves, but each of the 
four will interact with the bulk neutrino with its own coupling constant. In this model, the 
observed mixing between active neutrinos is induced from their individual couplings with 
the bulk. The fourth neutrino is forced on us by the strong coupling, because in that limit 
one brane neutrino always disappears into the KK towers. The details of how that comes 
about will be reviewed in the next section. In order to have three mass eigenstates on the 
brane in the presence of a strong coupling, corresponding to the three active neutrinos with 
definite masses, we must start from four brane neutrinos in the no-coupling limit. We shall 
refer to these four brane neutrinos before the coupling is turned on as 'flavor' neutrinos. 
Three of them are the usual active neutrinos, produced and absorbed by weak interactions 
in the usual way, but of course not the fourth because it is sterile. In other words, the 
fourth flavor is really flavorless. It is a flavor(-\ess) state only in the sense that a hadron is 
a co/or(-singlet) state. 

A distinct feature of the strong coupling is the absorption of the active neutrino fluxes by 
the bulk. This absorption is the result of a destructive interference of the infinite number of 
amplitudes oscillating through the KK mass eigenstates. If the coupling were weak, only a 
few nearby KK modes are reached. In that case more complicated oscillatory patterns (than 
those without the presence of any KK mode) are setup, but the amplitude remains oscillatory 
in nature, in the sense that it will return to its full flux at some appropriate time. When the 
coupling gets stronger, more KK modes are involved, the destructive interference between 
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these modes begins to take shape to dampen the oscillation amplitude over a finite time 
period. This feature can be seen in the numerical simulations in the paper of Barbieri et al 
and the paper of Cosme et al in Ref. [2]. Finally, as the coupling strength becomes infinite, all 
KK modes participate, destructive interference from these modes becomes complete, every 
hint of oscillation from these modes is wiped out, and the bulk becomes purely absorptive. 
Neutrino oscillations now proceed completely through the three mass eigenstates on the 
brane, but their amplitudes are now damped by the absorption from the bulk. 

There are nine possible oscillation channels, from any one of the three active neutrinos 
to any other one. In the MM, the absorptions from these nine channels are described 
by a single parameter m 4 . In other words, there are eight relations between the amount 
of absorptions in these nine channels. This distinct feature of absorption can be used to 
separate this mechanism of loss of amplitude from others, such as decay, as we shall discuss 
later. For now, let us describe in more detail what this parameter m 4 is, in what sense the 
MM is minimal, and what is meant by the strong coupling limit. This model starts with 8 
real parameters, the four Majorana masses m a of the flavor brane neutrinos, and their four 
Dirac-mass couplings d a to the neutrino in the bulk. There is a single bulk neutrino without 
a five-dimensional mass, whose KK energy spectrum consists of the integers, in the unit of 
1/R. The overall coupling strength is measured by d 2 = J2 a =i The strong coupling limit 
is taken to mean that this strength is much larger than the other dimensionless parameters 
involved, namely, d ^> m a , e a = d a /d, and 1. If the Dirac mass d/R is similar to the electron 
mass, then since the Majorana masses are small, we only require the Dirac mass to be much 
larger than the KK mass gap. This condition d 3> 1, with d/R ~ 0.5 MeV, is equivalent 
to R 3> 1/(0.5 MeV) ~ 400 fm, so it remains to be in the strong coupling limit even to 
fairly small extra dimensions. In that limit, we are left with 7 parameters. Six of them can 
be taken to be the active neutrino masses M 1 , M 2 , M 3 , and the mixing angles #i 2 , O23, and 
#i 3 . The 7th parameter m 4 , which is the Majorana mass of the sterile neutrino before any 
coupling is turned on, measures the amount of absorption by the bulk. Absorption is absent 
in the limit m 4 — > 00. 
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We now argue that an experimental observation of these distinct absorption patterns 
is a good indication of the presence of extra dimensions. As discussed above, absorption 
occurs only when many KK modes participate in a destructive interference. Since there is no 
simple reason to expect the presence of many sterile neutrinos without an extra dimension, 
we think that the presence of extra dimensions is strongly indicated if such an observation 
is made. However, partial loss of an oscillation amplitude may also come from decay [6], 
oscillation into a very heavy neutrino [7], or other mechanisms. To distinguish these other 
possibilities from the absorption by the bulk, a quantitative analysis is required. In the 
absorption mechanism, a single parameter 7714 governs the nine oscillation channels, so one 
can predict the amount of absorption of the other eight channels when one is measured. If 
this is verified, then one can be reasonably sure that the loss is due to the absorption in 
the bulk, and not something else. In other words, the presence of an extra dimension is 
indicated. 

In this paper we shall show that the results of the MM are rather generic. These results 
include the control of all absorptions by a single parameter m 4 , which allows us to distinguish 
extra dimensional absorption from others. They also include the natural connection between 
a small solar to atmospheric mass-gap ratio and the smallness of the reactor angle #13. We 
will show that these features are preserved in a much larger class of models in which the 
bulk spectrum is arbitrarily altered, so that the nature and the details in the bulk are not 
important, only its presence is. This is so because the amount of absorption in the strong 
coupling limit is controlled by the same seven parameters in the MM, and not the details of 
the KK spectrum. 

In Sec. 2 we discuss the mass matrix of the MM, and its generalization. The diagonaliza- 
tion of this infinite-dimensional matrix to find the mass eigenvalues and the PMNS mixing 
matrix is worked out. These results are used in Sec. 3 to compute the oscillation amplitudes 
and the absorption by the bulk. A conclusion can be found in Sec. 4. More complicated 
technical aspects are discussed in Appendix A. Throughout the rest of this article, we will 
not restrict the number of brane neutrinos to / = 4. 
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II. THE STRONG-COUPLING MODELS 



The mixing of four- dimensional neutrinos is given by a (real symmetric) mass matrix 
M. In the Minimal Model (MM), it is 

/ m D\ 



where m = diag(mi, rri2, • • • , rrif) is the Majorana mass matrix of the / brane neutrinos. 
B = diag(0, +1, — 1, +2, — 2, • • •) is the mass matrix of the bulk neutrinos, which are the 
Kaluza-Klein (KK) states of a single massless neutrino in a five-dimensional bulk. The 
coupling between the brane and the bulk neutrinos is assumed to be of the Dirac-mass type, 
given by the / xoo matrix D, with D an = d a . The matrix elements of D are independent of 
n G Z because every brane neutrino couples to the whole neutrino in the bulk, so it couples 
equally to each of its KK states. 

In this paper we generalize B to an arbitrary real symmetric matrix. In other words, 
we assume the spectrum of the bulk states that couple to the brane neutrinos to be rather 
arbitrary. They no longer need to come from a single massless neutrino in the bulk, nor do 
the neutrinos in the bulk need to be restricted to a five-dimensional spacetime. Some of the 
concrete models within this category will be discussed in Sec. 3. 

We shall use the discrete index n to label the bulk states. In the MM, n G Z, but in 
general, the set of n can be quite complicated. We shall denote it by K . 

The eigenvalue equation for the mass matrix is 



where w is a /-dimensional column vector with components w a , and v is an oo-dimensional 
column vector with components v n . A is the eigenvalue. In component form, (2) reads 




(1) 




(2) 



m a w a + d a A 



Xw, 



(3) 



b + (Bv) 



n 



Xv, 



n ■ 



(4) 



where 
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A = E u «> 

n 
f 

b =J2 d a W v ( 5 ) 
o=l 

We shall choose the normalization of the eigenvector so that 6=1. 

Since B is real symmetric, it can be diagonlized by a real orthogonal matrix O, so that 

B = 0~ 1 -fj,-0, (6) 

where fj, is a diagonal matrix with real matrix elements /i n . In the MM, O mn = 5 mn and 
H n = n. Defining u = Ov, and remembering the normalization 6 = 1, (4) becomes 

£ n + n n u n = Xu n , (7) 

where 

Cm = E Orrin- (8) 



n 



The eigenvector components can be solved from (3) and (4) to be 

A - ji n 



w a = A-^ = (Ad)-^- (9) 
A - m a X-m a 



For later convenience we have expressed d a = de a , where 



/ 

a=l 

i = E4 (io) 



a=l 

The constant A may now be computed to be 



A = J2 V n = E °mnU m = £ y^-, (H) 



n m.n 



where orthogonality of the matrix O has been used. 

The eigenvalue equation is obtained from (5) and (9) and the normalization condition 
6 = 1 to be 



1 = E d a w a = Ad 2 E t— ^ — = Ad " r - (12) 

„ 1 „ i A — Tin 

a=l a=l a 

In the strong-coupling limit when d — > oo, with the other parameters e a , m a , fi k and ^ kept 
fixed, eigenvalues either satisfy A(X) = 0, or r(A) = 0. We refer to the former as bulk 
eigenvalues, and the latter as brane eigenvalues. 

It is impossible to solve these eigenvalues analytically, especially for large /, but we 
know enough about them to make a substantial progress. Let us first consider the brane 
eigenvalues determined by r(A) = 0. Since r(A) goes to +oo at A = m a + and — oo at m a — , 
it must cross zero somewhere between one m a and the next. Consequently there are / — 1 
brane zeros all together, sandwiched between consecutive pairs of m a 's. We shall denote 
these brane eigenvalues as M i: with % = 1,2, •••,/ — 1. Note that these eigenvalues are 
completely independent of what B is. In particular, they are identical to those of the MM 
with the same parameters m a and d a . 

The bulk eigenvalues behave similarly because A(X) has the same structure as r(A). 
Hence there is one bulk eigenvalue between each consecutive pairs of /i n 's. The set of bulk 
eigenvalues will be denoted by K. In the MM, /i n = n G Z = K , and the bulk eigenvalues 
are K = Z + \. 

We proceed to consider the mixing matrix U. Its columns are the normalized eigen- 
vectors, with components U a \ = w a /N and U n \ = v n /N. The norm iV 2 of the original 
eigenvector (w a ,v n ) is given by iV 2 = (Ad) 2 s + T, where 

f P 2 

E ' 



a =i (A-m a ) 2 ' 

T W = E U n = E o kn u k o en u e 



n,k,e 



= ?(^' (13) 

If d is large but not infinite, the eigenvalues will shift somewhat, but they are still 
bounded between consecutive m a 's or consecutive /x n 's. The normalization factors in (13) 



have a d dependence, to be denoted by s<f(A) and T^(A). By definition, Soo(A) = s(A) and 
T oc (A)=T(A). 

III. TRANSITION AMPLITUDES IN THE ABSORPTIVE BULK 

Using the results obtained in the last section, we can calculate the transition amplitude 
A a b from a brane neutrino of flavor b and energy E, to a brane neutrino of flavor a after 
it has traversed a distance L = 2Et/U 2 . The transition amplitude is determined by the 
formula 

A ab (r) = £ KxU b xe- ix2r = A s ab (r) + A«(r), (14) 

A 

where A s is the contribution from the brane eigenvalues A = M i? and A K is the contribution 
from the bulk eigenvalues A G K. We shall see that A^ b (r) for r > is zero in the strong 
coupling limit, thus rendering the bulk an absorber of the brane neutrinos. This is the 
signature for the presence of an extra dimension. 

When d — > oo, the quantities v n ,w a /Ad, s and T are all of order 1, so the magnitude 
of w a is determined by Ad and the magnitude of N 2 is determined by (Ad) 2 . According to 
(12), Ad = l/(dr). For bulk eigenvalues, r = 0(1), so Ad = 0(l/d). This implies N 2 ~T 
and U a \ — 0(l/d). In that case the bulk components of an eigenvector are much larger than 
the brane components. For brane eigenvalues, A = 0(1), hence Ad = 0(d) and w a = 0(d). 
In that case the brane components of an eigenvector dominate and iV 2 ~ (Ad) 2 s. 

Let us denote the large-rf value of U a Mi by V ai , and l/(Aj — m a ) by x ai . Then 

V ai = e a x ai /yfc (1 < a < /, 1 < i < f - 1), (15) 

and 

A S ab (r) = EOfce-**?'. (16) 
i=i 

As it stands, V is a / x (/— 1) matrix, but we can make it into a square fxf matrix by letting 
the last column to be V a f = e a . The meaning of this last column will be discussed later. 
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Note that we can write 14/ in the same form as the other V ai , namely, V a f = e a x a f / y/sj, 
provided we let A = oo. The resulting square matrix 

e 2 x 2 i/^fs{ e 2 x 22 /\/s2~ ■■■ e 2 x 2 j-i/y/sf^ e 2 



V 



(17) 



\efx n ly/sl e/x/s/v^ ••• CfXfj-J ' y/sJH e f J 
is real orthogonal. This will be shown in Appendix A. 

Eq. (17) contains the parameters e a ,m a , and Mj, but they are not all independent. The 
/ e a 's are related by (10), so only / — 1 of them are independent. The brane eigenvalues 
Mj are solutions of r(A) = 0, so they are determined by the 2f — 1 parameters e a and m a . 
For large /, an analytic solution does not exist, so Mj cannot be expressed analytically 
in terms of e a and m a , which makes any computation involving Mj, such as (17) and (16), 
difficult to do in a closed form. To overcome this difficulty, it is crucial to make the following 
observation. 

Instead of the / — 1 parameters e a , we shall adopt the / — 1 values of Mj as independent 
parameters. Then e 2 a is a rational function of Mj and m a , as we shall show below. In this 
way everything can be carried out analytically, and this fact is used heavily in Appendix A 
to show the orthogonality of the matrix V. 

The crucial point is the observation that r(A) is a meromorphic function of A, with / 
simple poles at A = m a , and / — 1 zeros at A = Mj. Moreover, it follows from the definition 
of r in (12) and the constraint on e a in (10) that r(A) — > 1/A for large A. Hence from 
complex variable theory we can conclude that 

r(A) = Jf^A. m 

Since e 2 a is the residue of r(A) at the simple pole m a , it follows that 

2 _ UkZ\(m a - M k ) _ Uk=l(-Xak) nq x 

where for later convenience we introduce the abbreviations 
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x ai = Mi- m a , 
mab = m a - m b , 

Mij = Mi- Mj. (20) 

If there is one zero Mj between each consecutive pair of poles, then el > 0. We shall refer 
to this as the physical range. Otherwise, e a < may occur. 

Putting (19) into (13), we can compute Sj = s(Mj). It is shown in Appendix A that the 
sum over a can be carried out to yield the simple result 

* - -f^. (2D 

1 lc=l ^ci 

We consider now the contribution from the bulk eigenvalues. Since U a \ is unitary, it 
follows from (14) that A ab (Q) = 5 ab , hence 

A«(0) = S ab -A s ab (0). (22) 

Using (16) and the unitarity of the matrix V in (17), we conclude that 

Ai(0) = 5 ab - V a * f V bf = 5 ab - e a e f . (23) 

Therefore 

A« b (0) = e a e b . (24) 

The contribution from the bulk eigenvalues can also be obtained directly from (14) and 
the paragraph following that equation to be 

A K (r\ = V 1 ea6b r - aV (oft 

ah ^ K (dr)iT d (\)(\-m a )(\-m b ) e " ^ 

For a fixed A, both r and T are of order 1 as d — > oo, so the contribution from each bulk 
eigenvalue to the sum is 0{l/d 2 ). Since there is an infinite number of bulk eigenvalues, the 
total contribution to the sum in (25) is not necessarily zero. In fact, we know from (24) that 
Aa b (0) remains finite at infinite d. 
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There is no harm in dropping a finite number of terms in the infinite sum since each 
term contributes 0(l/d 2 ). For that reason we may assume A in the sum to be much larger 
than all the m a 's, in which case r(A) ~ 1/A 2 because of (10), so (A — m„)(A — m b )r 2 (\) ~ 1. 
The sum (25) can therefore be simplified to 

^ab( T ) = e a£b J/ U exp(-iAV) = e a e b g(r), (26) 

where 

9(r)=T,^7T-^M-^ 2 r). (27) 

Since g(0) = 1, we must have Z)agk^T 1 (^) — d 2 for large d. Coupling this with the earlier 
observation that T^(A) = 0(1) when d — > oo, we conclude that T d " 1 (A), as a function of A, 
is of order 1 until |A| = 0(d 2 ) ) after which it decreases sufficiently fast in A for the series to 
converge to d 2 . Assuming /j, n to be smooth, and the distances between one [i n and another 
are bounded from above, the additional oscillatory factor exp(— i\ 2 r) is going to render 
g{r) — > for finite r when d — > oo. This can be seen by using the natural variable y = \/d 2 . 
A very explicit calculation for the MM based on this reasoning can be found in Ref. [3] and 
[4]- 

In conclusion, the bulk contribution is 

A« b {r) = e a e b g(r), (28) 
with g(0) = 1 and g(r) ~ for r > 1/d 4 . 



IV. CONCLUSION 

In two previous publications [3,4], a five- dimensional minimal model (MM) of neutrino 
oscillation was discussed. It is a model with a strong brane-bulk coupling. It contains 
seven parameters, six of which can be taken to be the three masses of the active neutrinos 
and the three mixing angles. The seventh, WI4, controls the absorption by the bulk of 
the oscillation amplitudes in the nine channels. The MM model predicts the smallness of 
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the solar to atmospheric mass-gap ratio and the smallness of the reactor mixing angle are 
naturally related. We argued that this feature of the bulk absorption controlled by a single 
parameter 7714 can be used to detect an extra dimension. In this paper, we showed that both 
of these features depend on the presence of one or more extra dimensions, but not on the 
details of them, in the sense that they do not depend on the details of the KK spectrum. 
With this finding, we may now regard them to be generic features of strongly coupled extra 
dimensional models, and not just for the MM. 

This research is supported by the Natural Sciences and Engineering Research Council of 
Canada and by the Fonds de recherche sur la nature et les technologies of Quebec. 



APPENDIX A 

Several formulas quoted in the text will be proved in this appendix. 

To prove the formula (21) for = s(Mj), we start from its definition in (13): 

-=Sc«^ s S^- (A1) 

Considered as a function of Mj, Si has a pole at every m a . These are simple poles because 
e 2 a has a zero at the same point. The residue at Mj = m a can be obtained with the help of 
(19) to be — Il{^!(— x ak)/ n{^ a m ac- Moreover, because of (10), s$ asymptotically appraches 
1/Mf at large Mj. The function on the right hand side of (21) has exactly the same poles, 
the same residues, and the same asymptotic behavior, hence the difference between that 
function and Sj is an entire function which vanishes at infinity. This entire function must be 
zero, hence Sj is given by the function in (21). 

Next, we want to show the matrix fxf V in (17) to be an orthogonal matrix. We shall 
do that by showing each row of V to be normalized, and two different rows to be mutually 
orthogonal. 

Let us first prove that the ath row and the 6th row are orthogonal, when a 7^ b. The dot 
product of these two rows is 
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e a e b 



/-i 

E 



i 



= e a e b [1 + E\ 



r{ (Mi - m a )(Mi - m b )si 
For this to be zero we must have E — — 1. We will now show why this is so, again using the 
complex variable theory. 



Recall from (21) that 



ULi(Mi-m c ) 



Si = 



Hence 



E = _ , ^H^aj b (M i -m c ) 



i nfcW-M*) 

Consider E as a function of Mi. It has a simple pole at Mi = Mj for all j > 1. Such a pole is 
contained in the i = 1 term and the % — j term of the sum. The residue at Mi = Mj coming 
from the i — 1 term is — Y[{^ b {Mj — m c ) / n{^lj(Mj — M k ), and the residue coming from 
the i — j term is exactly the opposite. Hence the total residue is zero and E is an entire 
function of M^ Asymptotically at large M 1; it follows from (A2) that the only contribution 
comes from the % — 1 term, which gives E — — 1, hence E is the constant function — 1, which 
is required to show the vanishing of (A2) and the orthogonality of two different rows of the 
matrix V. 

Next, we show that every row of the matrix V is normalized, namely, 

1 



/-i 

i + E 



1 (Mi - m a ) 2 Si 



= 1. 



(A2) 



We will show this in the form 



i + E 



i 



^ (M - m a ) 2 Si 
Using (19), (21), this is equivalent to showing that 



o2 ' 



(A3) 



(A4) 



Let us label the function on the left as L(m a ), and the function on the right as R(m a ). 
The function L(m a ) goes to 1 for large m a , and so does R{m a ). Moreover, L{m a ) has a 
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simple pole at m a = A« for every i, with residue Y\{^ a X&/ Ylj^l M^. The function R(m a ) 
also has a simple pole at m a = Aj for every i, with the same residue. Hence L{m a ) — R(m a ) 
is an entire function that vanishes at infinity, so it is the constant function 0. This proves 
(A4), and hence the normalization condition (A2). 

We have thus shown that the matrix V is a real orthogonal matrix. 
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